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We develop an understanding of the anomalous metal state of the parent compounds of recently 
discovered iron based superconductors starting from a strong coupling viewpoint, including orbital 
degrees of freedom. On the basis of an intermediate-spin (5'=1) state for the Fe'^'*' ions, we derive a 
Kugel-Khomskii spin-orbital Hamiltonian for the active t2g orbitals. It turns out to be a highly com- 
plex model with frustrated spin and orbital interactions. We compute its classical phase diagrams 
and provide an understanding for the stability of the various phases by investigating its spin-only 
and orbital-only limits. The experimentally observed spin-stripe state is found to be stable over 
a wide regime of physical parameters and can be accompanied by three different types of orbital 
orders. Of these the orbital-ferro and orbital-stripe orders are particularly interesting since they 
break the in-plane lattice symmetry - a robust feature of the undoped compounds. We compute 
the magnetic excitation spectra for the effective spin Hamiltonian, observing a strong reduction of 
the ordered moment, and point out that the proposed orbital ordering pattern can be measured in 
resonant X-ray diffraction. 
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I. INTRODUCTION 

The begining of this year marked the discovery of a 
new and very unusual family of high temperature super- 
conductors: the iron pnictides. Superconductivity at 26 
K was discovered in fluorine doped rare-earth iron oxyp- 
nictide LaOFeAs^i^. In subsequent experimental studies 
involving different rare earth elements a superconduct- 
ing Tc larger than 50 K was reported^i^i^. Since then a 
large number of experimental and theoretical papers have 
been published, making evident the immense interest of 
the condensed matter community in this subject^. 

It has become clear that the iron pnictide superconduc- 
tors have, besides a number of substantial differences, 
at least one striking similarity with the copper oxides: 
the superconductivity emerges by doping an antiferro- 
magnetic, non-superconducting parent compound. This 
antiferromagetism is however of a very unusual kind. In- 
stead of the simple staggered '(tTjTt)' antiferromagnetism 
of the undoped cuprates, this 'stripe' or '(tt, 0)' spin or- 
der involves rows of parallel spins on the square Fe-ion 
lattice that are mutually staggered^. In fact, before this 
order sets in a structural phase transition occurs where 
the two in-plane lattice constants become inequivalent. 
This structural distortion is very small, but it appears 
that the electron system undergoes a major reorganiza- 
tion at this transition. This is manifested by resistivity 
anomalies, drastic changes in the Hall- and Seebeck co- 
efficients, and so forth^. Although the magnetic- and 
structural distortion appear to be coincident in the 122 
family^i^, in the 1111 compounds they are clearly sepa- 



rated , and there it is obvious that the large scale changes 
in the electron system occur at the structural transition 
while barely anything is seen at the magnetic transition. 

Given that the structural deformation is minute, this 
is an apparent paradox. Assuming that only the spins 
matter one could envisage that the spin ordering would 
lead to a drastic nesting type reorganization of the Fermi 
surfaces, causing a strong change in the electronic prop- 
erties. But why is so little happening at the magnetic 
transition? One could speculate that the spins are fluc- 
tuating in fanciful ways, and that these fluctuations re- 
act strongly to the structural chang o^^d^d^ . Such pos- 
sibilities cannot be excluded on theoretical grounds but 
whichever way one wants to proceed invoking only spins 
and itinerant carriers: one is facing a problem of princi- 
ple. 

This paper is dedicated to the cause that valuable 
lessons can be learned from the experiences with man- 
ganites when dealing with the pnictides. A crucial lesson 
learned over a decade ago, when dealing with the colos- 
sal magneto resistance (CMR) physics of the mangan- 
ites, was the demonstration by Millis, Littlewood, and 
Shraimau"'^'?' that the coupling between fluctuating spins 
and charge carriers can only cause relatively weak trans- 
port anomalies. In the pnictides one flnds that the re- 
sistivity drops by a couple of milliohm centimeters, that 
the Hall mobility increases by 2-3 orders of magnitude, 
and most signiflcantly the Seebeck coefflcient drops by 
an order of magnitude from a high temperature limit or- 
der value of 40 ji V/K in crossing the transition. It is 
very questionable if spin-carrier coupling of any kind, be 
it itinerant or strongly coupled, can explain such large 
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changes in the transport properties. 

A. Role of Electron-Electron Interactions 

Comparing the pnictidcs with tlie cuprate supercon- 
ductors there is now a consensus that in two regards 
these systems are clearly different: (i) in the pnictide 
system no Mott insulator has been identified indicating 
that they are 'less strongly correlated' than the cuprates 
in the sense of the Hubbard type local interactions; (ii) in 
the pnictide one has to account for the presence of several 
3d orbitals playing a role in the low energy physics, con- 
trasting with the single 3(ij.2_y2 orbital that is relevant 
in the cuprates. 

As a consequence, the prevailing viewpoint is to regard 
the pnictides as LDA metals, where the multi-orbital 
nature of the electronic structure gives rise to a multi- 
sheeted Fermi surface, while the 'correlation effects' are 
just perturbative corrections, causing moderate mass en- 
hancements and so on. 

Although there is evidence that the system eventually 
discovers this 'Fermi- liquid fixed point' at sufficiently low 
temperatures, it is hard to see how this can explain the 
properties of the metallic state at higher temperatures. 
The data alluded to in the above indicate pronounced 
'bad metal' behavior, and these bad metal characteris- 
tics do not disappear with doping. In fact, one can ar- 
gue that the term 'bad metal' actually refers to a state 
of ignorance: it implies that the electron system cannot 
possibly be a simple, coherent Fermi-liquid. 

B. Spin-Charge-Orbital Correlations 

Another important lesson from the manganites is that 
the presence of multiple orbitals can mean much more 
than just the presence of multiple LDA bands at the 
Fermi-energy. Manganite metals have a degree of itin- 
eracy in common with the pnictides, but they still ex- 
hibit correlated electron physics tied to orbital degen- 
eracy which is far beyond the reach of standard band 
structure theory. 

The seminal work by Kugel and Khomskii in the 1970's 
made clear that in Mott insulators orbital degrees of 
freedom turn into dynamical spin like entities that are 
capable of spin-like ordering phenomena under the con- 
dition that in the local limit one has a Jahn- Teller (or- 
bital) degeneracy^''. The resulting orbital degrees of free- 
dom can have in dynamical regards a 'life of their own'. 
This manifests itself typically in transitions characterized 
by small changes in the lattice accompanied by drastic 
changes in the electronic properties. 

In the manganites there are numerous vivid examples 
of the workings of orbital orderingi^^i^iii. Under the 
right circumstances one can find a transition from a high 
temperature cubic phase to a low temperature tetrago- 
nal phase accompanied by a quite moderate change in the 



lattice, but with a change in the electron system that is 
as drastic as a 'dimensional transmutation': this system 
changes from an isotropic 3D metal at high temperature 
to a quasi 2D electron system at low temperatures where 
the in-plane resistivity is orders of magntitude lower than 
the c-axis resistivit y "'^^i^^'^" . 

The explanation is that one is dealing in the cubic man- 
ganite with a Mn'^^ ion with an Cg Jahn- Teller degener- 
acy involving "idx^-y^ and Zd^^'^-i orbitals. In the low 
temperature 'A-phase' one finds a 'ferro' orbital order 
where cooperatively the x^-y^ orbitals are occupied. This 
greatly facilitates the hopping in the planes while for sim- 
ple orthogonality reasons coherent transport along the 
c-axis is blocked. Since the d-electrons only contribute 
modestly to the cohesive energy of the crystal, this large 
scale change in the low energy degrees of freedom of the 
electronic system refiect only barely in the properties of 
the lattice. On the other hand, this orbital order is a 
necessary condition for the spin system to order, and at 
a lower temperature one finds a transition to a simple 
staggered antiferromagnet, in tune with the observation 
that in the A-phase the effective microscopic electronic 
structure is quite similar to the ones found in cuprate 
planes. 

The ruthenates are another class of materials in which 
the orbital degrees of freedom play a decisive role, in both 
the metallic and insulating phases. Bilayer Ca3Ru207, 
for instance, has attracted considerable interest because 
the observed CMR-effect is possibly driven by orbital 
scattering processes among the conduction electrons^^'^^. 
Another example is TI2RU2O7, in which below 120 K its 
3D metallic state shows a dramatic dimensional reduction 
and freezes into a quasi- ID spin system, accompanied by 
a fundamental orbital reorganization^^iSi. 

It is very remarkable that the groundstate of all 
iron pnictides is characterized by a very similar spatial 
anisotropy of the magnetic exchange interactions: along 
one direction in the plane the Fe-Fe bonds are strong and 
antiferromagnetic, whereas in the orthogonal direction 
they are very weak and possibly even ferromagnetic^^. 
With all the others, also this observation is consistent 
with our hypothesis that the undoped iron pnictides are 
controlled by 'spin-charge-orbital' physics, very similar 
in spirit to the ruthenates and manganites. 



C. Organization of this Paper 

In Sec. II of this paper we derive the spin-orbital 
Hamiltonian starting with a three-orbital Hubbard model 
for the iron square lattice of the iron pnictides. The phase 
diagrams in the classical limit of this Hamiltonian are 
discussed in Sec. HI. We analyze the various phase tran- 
sitions by also considering the corresponding spin-only 
and orbital-only models. Sec. IV deals with the results 
on magnetic excitation spectra, which provide a possi- 
ble explanation for the reduction of magnetic moment, a 
central puzzle in the iron superconductors. We conclude 
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by commenting briefly on how the itineracy may go hand 
in hand with the orbital 'tweed' order that we put for- 
ward in the present study, and point out that the 'tweed' 
orbital ordered state can, in principle, be observed in 
resonant X-ray diffraction experiments. 



II. SPIN-ORBITAL MODEL FOR IRON PLANES 
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As stated above, the superconducting iron-pnictides 
are not strongly coupled doped Mott insulators. Staying 
within the realm of Hubbard-model language they are 
likely to be in the intermediate coupling regime where 
the Hubbard J7's are of order of the bandwidth. To at 
least develop qualitative insight in the underlying physics 
it is usually a good idea to approach this regime from 
strong coupling for the simple reason that more is going 
on in strong coupling than in the weak coupling band 
structure limit. As the experience with for instance the 
manganites and ruthenates shows, this is even more true 
when we are dealing with the physics associated with or- 
bital degeneracy. The orbital ordering phenomena that 
we have already alluded to, take place in itinerant sys- 
tems but their logic is quite comprehensible starting from 
the strongly coupled side. 

Thus as a first step we will derive the spin-orbital 
model of pnictides starting from a localized electron 
framework. A condition for orbital phenomena to oc- 
cur is then that the crystal fields conspire to stabilize 
an intermediate spin {S = 1) ionic states. These crystal 
fields come in two natural varieties: one associated with 
the tetrahedral coordination of Fe by the As atoms, and 
a tetragonal field associated with the fact that the overall 
crystal structure consists of layers. When these crystal 
fields would be both very large the Fe 3c?^ ions would 
form a low spin singlet state. This is excluded by the 
observation of magnetism, and moreover band structure 
calculations indicate that the crystal fields are relatively 
small. 

The other extreme would be the total domination of 
Hund's rule couplings and this would result in a high spin 
S = 2 state, which appears to be the outcome of spin po- 
larized LDA and LDA-I-U calculations^. However, given 
that for elementary chemistry reasons one expects that 
the tetrahedral splitting is much larger than the tetrago- 
nal splitting there is the possibility that the Hunds rule 
overwhelms the latter but looses from the former, result- 
ing in an 'intermediate' 5 = 1 state. Although the issue 
is difficult to decide on microscopic grounds, for orbital 
physics to be relevant we need an intermediate spin state 
as in the present crystal field scheme this is the only ionic 
(f' state that exhibits a Jahn- Teller groundstate degen- 
eracy (see Fig. 1). 

In this situation the starting Hubbard model involves 
a non-degenerate \xy) and two doubly-degenerate \xz) 
and \yz) orbitals, as will be defined in subsection A. The 
details of the derivation of the model are given in sub- 
section B. The derivation does not assume any specific 
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FIG. 1: (Color online) (a) Fe square lattice (black circles) 
and relative positions of the As ions. The latter are located 
in adjacent layers above (filled red squares) and below (empty 
red squares) the Fe plaquettes. (b) Schematic illustration of a 
ground-state configuration of the Fe ions corresponding to 
an intermediate S — 1 spin state, (c) Multiplet structure of 
the d| djd| charge excitations for localized Cg electrons. 



structure for the hopping parameters and hence, is com- 
pletely general. The algebra involved in the derivation 
is tedious but straightforward and a general reader may 
wish to skip subsection B and jump directly to subsec- 
tion C where we discuss the relevant hopping processes 
for the Fe-As plane. Incorporating these hopping param- 
eters leads to the model relevant to the iron plane. 



A. Hubbard model for pnictide planes for the 
intermediate-spin state 

The iron ions are in a dP configuration where we assume 
the low lying Cg orbitals to be fully occupied due to a 
large crystal-field splitting between the eg and t2g states. 
The two remaining electrons occupy the three t2g orbitals 
\a) := \xz), \b) := \yz), and |c) := \xy) with x and y 
pointing along the bonds of the iron square lattice. Due 
to the Hund's coupling Jh between the t2g electrons, 
such a configuration leads to an 5 = 1 intermediate spin 
state of the Fe ions. Further, we incorporate a small 
tetragonal splitting A between the \xy) state and the 
\xz), \yz) doublet (see Fig. [1]). 

Assuming the eg electrons to be localized, the phys- 
ical situation is very similar to almost cubic vanadates 
like YVO3 or LaVOa where the two d-electrons of the 
V'^"'" ions occupy nearly degenerate t2g orbitals. Interest- 
ingly, in theses systems orbital ordering in the presence 
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of a small crystal-field splitting A can lead to C-type 
antiferromagnetism22i22ii2S characterized by an ordering 
wavevector Q = (tt, 7r,0). The effective Hubbard model 
for the t2g electrons consists of a kinetic energy part Tit , 
a crystal field sphtting Hd, and of the on-site electron- 
electron interactions Tiint, 



H = Tit + Tic 



Hi, 



(1) 



with a kinetic energy contribution that is much richer 
than in the vanadates. For the nearest neighbor bonds 
the effective hoppings between the Fe t2g orbitals have 
contributions from both direct d — d and d — p — d pro- 
cesses via As p-orbitals. These As ions are located in 
adjacent layers above or below the Fe ion plaquettes as il- 
lustrated in Fig.[T^. Because of this particular geometry, 
the indirect As mediated hoppings should be of similar 
strength for nearest and next-nearest neighbor Fe ions. 
At this point, we do not specify the effective hopping ma- 
trix elements t^^'-p between orbitals a, P = a,b, c along a 
particular bond and write the kinetic energy oper- 
ator in the most general form, 



(2) 



where d}^^ {diaa) creates (annihilates) an electron on site 
i in orbital a with spin a —'\, |. The crystal-field splitting 
between the t2g orbitals is simply given by 



with fiia = J2a '^laa and fli, 



7t 



In our case 



the electron energies are given by Cc = for the xy and 
= Cfc = A for the xz and yz orbitals. The electron- 
electron interactions are described by the on-site terms,— 



2 V 2 



with the Coulomb element U and a Hund's exchange el- 
ement Jh- 



B. Superexchange model 

In the limit of strong Coulomb repulsion, t ^ U, 
charge fluctuations d^dj ^ djd^ are suppressed and on 
each site the two t2g electrons have to form a state be- 
longing to the ground-state manifold of Tiint + Ticf in the 



two-electron sector. For sufficiently small crystal-field 
splitting, A^ < 8Jjj, these states are given by two 5=1 
triplets in which on each site either the xz or yz is unoc- 
cupied. This orbital degree of freedom can be viewed as 
a T = i pseudospin. From Eqs. ([3]), (|4]) we easily obtain 
Eq = U — 3Jh + a as the ground-state energy of the t^g 
sector. 

A general spin-orbital superexchange model can be de- 
rived by second order perturbation theory controlled by 
the kinetic energy contribution Tit, where we have to con- 
sider all virtual processes t^gt^g 



4-1 f3 



tigtig acting 



on the S — 1, T — 1/2 ground-state manifold. The most 
general superexchange Hamiltonian in the sense of Kugel 
and Khomskii for a given bond takes the form 



n 



KK 



EE4: 



^Bg^g. {Si, Sj 



(5) 



where S and T denote S — 1 spin and T — ^ pseudospin 
operators. The functional form of B only depends on 
total spins Si and Sj on the two sites in the intermedi- 
ate t\gt2g states. Whereas the single occupied site has 
necessarily s — 1/2 the other site can be in a high-spin 
(s = 3/2) or low-spin (s = 1/2) state. Likewise, the func- 
tions are determined by the pseudospins r^, Tj of 
the involved intermediate states. 

To derive the effective spin-orbital superexchange 
model we first have to find the multiplet structure of the 
virtual intermediate configurations. It is straightfor- 
ward to diagonalize Tid + 'Hint pi4p in the three-particle 
sector. The lowest energy we find for the quar- 



ts) tet of s ~ 3/2 high-spin intermediate states \ A2, |,s^ 
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energy is £(^^2) = E{'^A2) - 2Eq ^ U - 3Jh, where 
Eq = U — 3Jh -I- a is the groundstate energy in the 
sector. In order for the approach to be valid we have 
to assume that the system has a charge-transfer gap, 
U — 3Jh > and that the hopping matrix elements 
are sufficiently small compared to the charge-transfer 
gap. All the other multiplets consist of intermediate 
s = 1/2 doublets. The '^E multiplet with excitation en- 
ergy e(^E) — U consist of the two spin-i doublets 

\'E,^,a), = ^(2rfLrfL<-.-4X-.<a 



^E,-,a)2 



-di^Aadim 



dLd 



(6) 

rdim- (7) 



Finally, we have multiplets : which consist 

of spin-i doublets and invoke doubly occupied orbitals. 
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with excitation energies e(^Ti) = U and 6(^X2) = U+2Jh 
and 

\'T,%,l,<7h = dU^f^dl/^,Tv^dl/J\0) (9) 
\'T,%,\,<t)2 = 4.(^1 -4<T<iT'^T«)|0> (10) 



with Vzfz = Jh / y Jfj + ± + J|f)^ and excitation 
energies e(2T/)2) = [/ + T ^f^^TJ^. 

The resuhing charge-excitation spectrum is shown 
schematically in Fig. [TJ:. Although the single occupied 
t\g site of a virtual t\gt\g intermediate state gives no 
contribution to the on-site electron-electron interaction 
it can lead to an additional crystal-field energy A if the 
electron is in the a or & orbital. 

Let us first focus on the purely magnetic parts 
Bsi^s {Si, Sj) of the superexchange Hamiltonian, which 
can be determined entirely by group theoretical meth- 
ods. To be precise, we consider a two-ion system in the 
state 15*^, Ma) ® \Sb,Mb) which can be classified by the 
total spin St and the z-component A'lf Applying a hop- 
ping operator of the form Tit — —^^^{c^Aa^Bcr + h-c), 
which preserves the quantum numbers iS** and Mt be- 
cause of the spin-rotation invariance we obtain an inter- 
mediate state |sa,wa) ® |s_B,mB) with sa = Sa^ 1/2 
and SB = Sb i 1/2. The effective superexchange involv- 
ing intermediate spins sa, sb is given by the second order 
process 



E{Sf,SA,SB) = 



E 



\{sAmA, SBmBlHtlSiMt 
AE 



Using Clebsch-Gordan coefficients C^f'^^m — 
0ij2'7iim2 1 jm), we can express the total spin states as 

\StMt)^ E C^f^sMjSA,MA)<E>\SB,MB). 
Ma, Mb 

Since the operators and Ca are irreducible tensor oper- 
ators of rank 1/2 we can use the Wigner-Eckart theorem 
to obtain 



{sAmA\c\JSAMA) - \\c\\\C;;^ 
{sBmB\cB.\SBMB) - l|csrf;Jl",)„J-l)^-, 

where we have used || • || as a short-hand notation for 
the reduced matrix elements. Using these expressions 
we can rewrite the exchange energy as E{St, sa, sb) ~ 

'Ke^W'^^aW ' \\'^B\\)'^B{St,SA,SB), where we can express the 
function B in terms of a Wigner 6j-symbol as 



Sb{Sb + 1) + 2SaSb can be simplified further to 

^ -i2SA + l)i2SB + l)''{^'^^l^^'''^'2^ 



sign[(sA - S'a)(ss - Sb)]SaSb 



This expression we can evaluate for Sa ^ Sb = S — 1 
for the s = 3/2 high- and s = 1/2 low-spin intermediate 
states to obtain the (normalized) spin-projection opera- 
tors 



i?|,i(^.,5,) = -i(5,5, + 2) 



(11) 
(12) 



in agreement with Refs. l27l|28l . Hence, the Kugel- 
Komskii superexchange Hamiltonian for a given bond 
can be written as. 



-(5.S,-1)Q(2)(T.,T, 



(13) 



where Q^"-* are functions of orbital pseudospin operators. 
Their functional form can be obtained by tracking the 
orbital occupancies in the initial and final states during 
a virtual hopping process. In terms of spinless Fermi 
operators af, bf increasing the occupancy of the a or 
b orbital on site i the pseudospin-1/2 operators acting 
on the ground-state manifold can be expressed as T/ — 
{fiia - hib)/2, T+ = bfai, and T~ = afh, where fna = 
afai and fiH) = b^bi with the constraint Uia + fin, — 
1. Whereas it is straightforward to see that the general 
functional form is given by 



{f: + f/) + /(") if[ + ff) + 



(n) 



(14) 



it is quite tedious to determine the coefficients by acting 
with the hopping operator Tit (O on all states in the 
ground state sector and calculating the overlap of the 
resulting states projected on the different intermediate 
states listed above. The resulting explicit expressions 
are given in Appendix A. 



C. Hopping and Resulting Hamiltonian 



B{SuSa,sb) = -{2sa + 1){2sb + 1)\^^ \ 

SB OB Jt 



which by using the relation St{St + 1) = Sa{Sa + 1) + 



In the previous section we have derived the general KK 
superexchange Hamiltonian only assuming the effective 
hopping matrices to be symmetric, t^is = tpa. In or- 
der to write down the spin-orbital model specific to the 
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FIG. 2: Illustration of the effective hopping parameters tap, 
(a) between the dxz and dyz orbitals, and (b) those involving 
the dxy orbitals. The projections of the d^z and dyz orbitals on 
the Fe-plane are depicted in white and light grey, respectively, 
and the dxy orbitals are shown in dark grey 

pnictide planes we have to use the corresponding hop- 
ping parameters. We use the Slater-Koster integrals^ 
along with the geometry of the Fe-As planes to deter- 
mine all the hopping parameters involving the three 
orbitals on the nearest neighbor and next nearest neigh- 
bor Fe sites. This considerably reduces the number of 
independent hopping parameters that enter the Hamil- 
tonian. The direct d — d hoppings are considered to be 
much smaller therefore we use hoppings via the As-p or- 
bitals only which are given in Appendix B and depend 
on the direction cosines /, m, n of the As-Fe bond, as well 
as on the ratio 7 = (^pdn) / (pda)^'^^^. These result- 
ing effective hopping matrix elements between the t2g Fe 
orbitals are shown schematically in Fig. [2] and can be 
parametrized by the lattice parameter A = \n/l\ and 7. 

In Fig. [3] the dependence of the hopping matrix el- 
ements on the ratio 7 — (pdir) / (pda) is shown for a 
lattice parameter A = 0.7 which is slightly below the 
value resulting from the Fe-Fe spacing and the distance 
of the As ions to the Fe planes. Over a realistic range 
—0.2 < 7 < 0.2 we find a very strong dependence of the 
hopping amplitudes on 7 and therefore expect the sta- 
bility of possible phases to depend crucially on 7. This 
parameter cannot be obtained by geometrical consider- 
ations but depends for instance on how strongly the or- 
bitals delocalizc. 

Having specified the effective hopping parameters 
cki '■— ti/t between the Fe orbitals for nearest and next- 
nearest neighbors (see Fig. [2|) which are parametrized 
entirely by the ratio 7 = (jidn) / (pda) and the lattice pa- 
rameter A = \n/l\ we can now write down the effective 
KK model for the Fe planes. For convenience, we rewrite 
the Hamiltonian in the form 

nKK = Jj2 (U^^^^ + + (15) 

introducing an overall energy scale J — At^ /U . The or- 
bital bond operators are defined as^l — ^(Q'^' — Q'^^^) 
and f = - y|^(Q(^) -l-2g(2)) and depend on the effective 
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FIG. 3: (Color online) Various hopping parameters Qi ~ 
ti/t as illustrated in Fig. [2] as a function of the ratio 7 = 
(pd-K) / {pda) for the lattice parameter A = 0.7. 

couplings Ui , the relative strength of the Hunds coupling 
rj = Jh/U and the crystal-field splitting 6 = A/U. For 
the nearest neighbor bonds along x and y along the x±y 
diagonals the operators are given in Appendix C. 

III. CLASSICAL PHASE DIAGRAMS 

In this section we discuss the phase diagrams of the 
spin-orbital Hamiltonian in the classical limit. We have 
four parameters that enter the model: A and 7 deter- 
mine the relative strength of various hopping parameters 
and r] and S enter via the energy denominators. Zero 
temperature phase transitions are discussed in subsec- 
tion A, subsection B is devoted to the understanding of 
finite temperature transitions and subsection C analyses 
the phases in terms of the corresponding spin-only and 
orbital-only models. 

The results that we discuss below demonstrate that 
the Hamiltonian is highly frustrated in the spin and or- 
bital variables. While the spin frustration is largely due 
to the competing interactions between nearest- and next- 
nearest neighbors, the frustration in orbital sector is more 
intrinsic and exists within a single bond in the Hamilto- 
nian. The spin (tt, 0) state is found to be stable over 
a wide range of parameter space due to the strong nnn 
AF coupling. However, depending on the parameters, 
there are three possible orderings of the orbitals that ac- 
company the spin 'stripe' order. Two out of these three 
orbital ordering patterns break the in-plane symmetry of 
the lattice and hence are likely candidates for explain- 
ing the orthorhombic transition observed in the parent 
compounds. 

A. Zero temperature 

Since the effective KK Hamiltonian derived in section 
|TT]contains a large number of competing terms it is almost 
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FIG. 4: (Color online) 77-7 phase diagram for A = 0.7 and 
S = 0.01. ri = Jh/U, 7 = (pdn) / (pda) . The phases are 
denoted by their ordering wavevectors in the spin and orbital 
variables. Tz or Tx refers to the component of the orbital 
pseudospin that is saturated in the ordered state. 

impossible to anticipate what kind of spin-orbital order- 
ings are realized for different parameter values, in partic- 
ular since the signs and relative strengths of the effective 
hoppings ai between nearest and next-nearest neighbor 
Fe orbitals crucially depend on the ratio 7 = (pdi:) / (pda) 
as pictured in Fig. [3l While the parameters ai, and 
a-j do not show large relative changes over the range of 
7 shown in the figure, there are very clear crossings be- 
tween a2 and and as and ag. 

Recall that as and ag are the hoppings between nearest 
and next-nearest neighbors involving orbital |c) := \xy). 
If we infer the spin order arising purely from the non- 
degenerate |c) orbital, it suggests that the spin state 
should be (tt, 0)-ordered for ag < 2a| and (tt, 7r)-ordered 
otherwise. Therefore, this would imply that as 7 —0.2 
the magnetic superexchange resulting from the |c) or- 
bitals only favors (tt, tt) antiferromagnetism, whereas the 
(tt, 0) stripe AF becomes favorable for 7 0.2. 

A similar spin-only analysis for the degenerate orbitals 
I a), 1 6) is not possible and one has to treat the full spin- 
orbital Hamiltonian in order to find the groundstates. 
Nevertheless, the complicated variations in the hopping 
parameters already suggest that we can expect a very rich 
and complex phase diagram for the groundstate of the 
spin-orbital Hamiltonian. In particular in the region of 
intermediate 7 where the magnetic superexchange model 
resulting form the |c) orbitals only becomes highly frus- 
trated we expect the magnetic ordering to depend cru- 
cially on the orbital degrees of freedom. 

We first look at the classical groundstates of this 
model. We make use of classical Monte-Carlo method 
in order to anneal the spin and orbital variables simul- 



FIG. 5: (Color online) Schematic pictures of the three 
ground-state orbital ordering patterns that accompany the 
spin-stripe phase, (a) Orbital-ferro, (b) orbital-stripe, and 
(c) orbital-antiferro 



taneously, starting with a completely random high tem- 
perature configuration. Using this method we identify 
the various groundstates that exist for a combination 
of model parameters. In order to obtain a groundstate 
phase diagram, we minimize the total energy for a set of 
variational states which also include all the Monte-Carlo 
groundstates obtained for different choice of parameters. 

Fig. [4] shows the resulting T = phase diagram for 
varying 77 = Jh/U and 7 = [pd-K) / {pda) . The lattice pa- 
rameter A is fixed to 0.7, which is close to the experimen- 
tal value for the oxypnictides. The crystal-field splitting 
between the |c) and the |a),|6) orbitals is considered to 
be very small, 5 = A/t/ — 0.01. As expected, a large 
number of phases are present in the phase diagram. 

With increasing 7 we indeed find a transition from a 
(tt, tt) to a (tt, 0) antiferromagnet as suggested from the 
analysis of the frustrated magnetic superexchange model 
involving only the |c) orbitals. This is not surprising since 
the corresponding couplings a\ and/or a§ are sufficiently 
strong and as 7 0.2 the biggest hopping element is 
in fact given by ag between next-nearest neighbor |c) 
orbitals (see Fig. [3]). Whereas the (tt, 0) stripe magnet for 
large 7 is accompanied by an antiferro-orbital ordering 
of the Tz components corresponding to a checkerboard 
arrangement of the \a) and |6) orbitals (see Fig. [SJ;) for 
intermediate, small 7 we find two (tt, 0) magnetic phases 
possessing orbital orderings which are likely to break the 
inplane symmetry of the lattice structure. 

For small 77 we find a ferro-orbital arrangement of the 
Tz components corresponding to the formation of chains 
along the ferromagnetically coupled spin directions (see 
Fig. [5^). The existence of this orbital order crucially de- 
pends on the pre-existence of a spin stripe state, which 
generates magnetic-field-like terms for the orbital pseu- 
dospins. This will be discussed in detail when we try to 
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FIG. 6: (Color online) 77-7 phase diagram for A — 0.8 and 
S = 0.01. Note that the orbital ordered states that break the 
orthorhombic symmetry do not exist for this choice of A. 




0.00 0.05 0.10 0.15 0.20 0.25 0.30 
V 

FIG. 7: (Color online) rj-S/rj phase diagram for A = 0.7 and 
7 = {pdn) / {pda) = —0.05. This phase diagram illustrates the 
point that 5 is not a crucial parameter in the Hamiltonian. 



understand the thermal phase transitions. For larger rj 
the orbital order changes to an orbital-(7r, 0) 'tweed' pat- 
tern with a condensation of the components. This cor- 
responds to the formation of orbital zig-zag chains along 
the antiferromagnetically coupled spin direction as pic- 
tured in Fig. [5)3. Interestingly, the stripes in the magnetic 
and orbital sectors have the same orientation, contrary 
to the conventional Goodenough-Kanamori rules. How- 
ever, since we are dealing with a highly frustrated spin- 
orbital model involving nearest and next-nearest neigh- 
bor bonds these naive rules are not expected to hold. The 
'tweed' orbital order is expected to lead to a displace- 
ment pattern of the As ions, which can in principle be 
observed in X-ray differaction experiments. The 'tweed' 
orbital pattern should show up as a higher order struc- 
tural Bragg peak at (tt, 0). The orbital order might also 
be directly visible resonant X-ray diffraction at the iron 
K-edge, a technique that was pioneered in the mangan- 
ites^i^i^, and is nowadays available for all transition 
metal K-edges, in particular the iron one^^. Polariza- 
tion analysis and azimuthal angle dependence can distin- 
guish between charge, spin and orbital contributions to 
the resonant signal^^ which gives the possibility in the 
iron pnictides to single out the 'tweed' orbital pattern. 

The orbital stripe order persists to the regime of larger 
negative 7 where the magnetic order changes to the (tt, tt) 
antiferromagnet. This shows that the orbital 'tweed' 
state does not have spin-(7r, 0) order as a pre-requisite 
and therefore this orbital order can, in principle, exist 
at temperatures higher than the spin transition tempera- 
tures. In the regime of large Hund's coupling, r/ > 0.3 the 
system becomes ferromagnetic. This tendency is easy to 
understand since in the limit r/ — > 1/3 the charge-transfer 



gap closes and the KK model is dominated by processes 
involving the low-lying high-spin multiplet favoring 
a ferromagnetic superexchange. 

Let us further explore how the groundstate phase di- 
agram changes as we vary the lattice parameter A and 
the crystal-field splitting S. Fig. [6] shows the same phase 
diagram as in Fig.[4]but for a slightly larger separation of 
the As ions to the Fe-planes, A = 0.8. The two interesting 
phases with magnetic (tt, 0) and orbital-stripe and orbital 
ferro orderings do not appear in this phase diagram indi- 
cating that the stability of these phases crucially depends 
on the relative strength of nearest and next-nereast hop- 
pings which can be tuned by A. Presence of a tetracritical 
point is an interesting feature in this phase diagram. 

Finally, we analyze the dependence on the crystal field 
splitting S = A/U which so far we assumed to be tiny. 
Wc do not find any qualitative change of the groundstate 
phase diagram with increasing S. In particular, there 
are no new phases that appear and therefore the crystal- 
field splitting does not seem to be a crucial parameter. 
For example, the phase diagram in the ry-^/?7-plane for 
A — 0.7 and 7 — —0.05 shown in Fig. [7] indicates that 
a change in 6 only leads to a small shift of the phase 
boundaries. 



B. Finite temperature 

To obtain the transition temperatures for the various 
phase transitions, we track different order parameters as 
a function of temperature during Monte-Carlo annealing 
where we measure the temperature in units of the energy 
scale J. For example, the spin structure factor is defined 
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7V2 



(16) 



where {■■■)av denotes thermal averaging and N is the to- 
tal number of lattice sites. The orbital structure factor 
0(q) is defined analogously by replacing the spin vari- 
ables by the orbital variables in the above expression. 
Depending on the groundstate, different components of 
these structure factors show a characteristic rise upon 
reducing temperature. 

We fix (5 = 0.01, A = 0.7 and 7 = -0.05 and track the 
temperature dependence of the system for varying rj. For 
T = this choice of parameters corresponds to a cut of 
the phase diagram shown in Fig. |4] through four different 
phases including the two (tt, 0) stripe AFs with orbital 
orderings breaking the in-plane lattice symmetry. 

In Fig. [5] the temperature dependence of the corre- 
sponding structure factors is shown for representative 
values of the Hund's rule coupling 77. For small values 
of rj the groundstate corresponds to the orbital-ferro and 
spin-stripe state as shown in the phase diagram in Fig. 31 
Fig. [S^ shows the temperature dependence of ^(Tr, 0) and 
0(0, 0) which are the order parameters for the spin-stripe 
and orbital-ferro state, respectively. While the S'(7r, 0) 
leads to a characteristic curve with the steepest rise at 
T ~ 0.5, the rise in 0(0,0) is qualitatively different. In 
fact there is no transition at any finite T in the orbital 
sector. We can still mark a temperature below which a 
significant orbital-ferro ordering is present. The origin of 
this behavior lies in the presence of a Zeeman-like term 
for the orbital pseudospin. 

For 77 — 0.15, the phase diagram of Fig. |4] suggests a 
state with stripe ordering in both spin and orbital vari- 
ables. We show the temperature dependence of S'(7r,0) 
and O{tt,0) in Fig. [8)3. In this case both the spin and 
orbital variables show a spontaneous ordering, with the 
spins ordering at a much higher temperature. An inter- 
esting sequence of transitions is observed upon reducing 
temperature for r] = 0.18 (see Fig. |8}:). This point hes 
close to the phase boundary between spin-stripe and spin- 
ferro state with the orbital-stripe ordering. The spin- 
stripe order parameter ^(Tr, 0) shows a strong rise near 
T = 0.4. The orbital stripe order sets in at T ~ 0.15. 
The onset of this orbital order kills the spin-stripe or- 
der. Instead, we find that the (tt, tt) components of the 
spin structure factor shows a strong rise along with the 
(tt, 0) component of the orbital structure factor. Finally 
for r/ — 0.24 the orbital stripe ordering is accompanied 
by the spin antifcrro ordering, with the orbital ordering 
setting in at slightly higher temperatures (see Fig. [SJi) . 

The results shown in Fig. [S] are summarized in the 
T — rj phase diagram shown in Fig. [9l For small 77, 
the groundstate is spin-stripe and orbital-ferro ordered. 
While the spin order occurs at higher temperatures, there 
is no genuine transition to the orbital-ferro state. The 
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FIG. 8: (Color online) Relevant structure factors as a function 
of temperature for different values of ij. The lattice parameter 
and the relative strength of a and n hopping are fixed as 
A = 0.7 and 7 = —0.05, respectively. 



orbital-ferro state is driven by the presence of a magnetic- 
field-like term for the orbital pseudospin in the Kugel- 
Khomskii Hamiltonian. The stability of the orbital-ferro 
state crucially depends on the presence of the spin-stripe 
order. The dotted line joining the black circles in the 
small-77 range is only to indicate the temperature below 
which the orbital-ferro order is significant. This typical 
temperature scale reduces with increasing 77, until the 
system finds a different groundstate for the orbital vari- 
ables. Note that the temperature scales involved are very 
small owing to the highly frustrated nature of the orbital 
model, nevertheless there is no zero-temperature transi- 
tion in this purely classical limit. 



The spin-stripe state remains stable with the transition 
temperature reducing slightly. The transition tempera- 
ture for the orbital stripe state increases upon further 
increasing 77. For 0.15 < 7; < 0.2, multiple thermal tran- 
sitions are found for the magnetic state. The spin-stripe 
order which sets in nicely at T ~ 0.35 is spoiled by the 
onset of orbital-stripe state, which instead stabilizes the 
spin (7r,7r) state. Beyond rj — 0.2, The orbital stripe 
state occurs together with the spin antiferro state, with 
the spin ordering temperatures slightly lower than those 
for the orbital ordering. For 7/ > 0.3, the spin state be- 
comes ferromagnetic. 
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FIG. 9: (Color online) T-rj phase diagram for S = 0.01, A 
0.7, and 7 = -0.05 



Corresponding Orbital-only and Spin-only 
models 



In an attempt to provide a clear understanding of the 
spin and orbital ordered phases, we derive the orbital 
(spin) model that emerges by freezing the spin (orbital) 
states. For fixed spin correlations, the orbital model can 
be written as 



+ E E T^tT^ + E (17) 

Here and below (•,•) and ((•,•)) denote bonds between 
nearest and next-nearest neighbor pseudospins on the 
square lattice, respectively, /i denotes the component 
of the orbital pseudospin. The effective exchange cou- 
plings for this orbital-only model are shown in Fig. (TU] 
as a function of the Hund's coupling 77 — Jh/U with 
the other parameters fixed as ^ = 0.01, 7 = —0.05, and 
A = 0.7, as before. The solid hnes are obtained by fixing 
the spin degrees of freedom by the classical ground-state 
configurations of the corresponding phases. For compari- 
son, the effective couplings for disordered spins are shown 
by dashed lines. 

Similarly, we can freeze the orbital degrees of freedom 
to obtain an effective Heisenberg model for spins, 

Ti-s = Jx ^ StSj + Jy ^ SiSj + Jd''^ SiSj. (18) 

The coupling constants Jx, Jy, and Jd for spins are plot- 
ted in Fig. [nj 
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FIG. 10: (Color online) The coupling constants as a function 
of 7; = Jh/U for the orbital-only model with frozen spin cor- 
relations for 5 — 0.01, 7 = —0.05, and A — 0.7. The couplings 
along X, y and diagonal directions are plotted in panels (a), 
(b) and (c) respectively. The single site term is plotted in (b) 
to indicate that this term arises due to a ferromagnetic bond 
along y-direction. The solid lines correspond to the ground 
state spin order and the dashed lines are for a paramagnetic 
spin state. The vertical dashed line indicates the location in rj 
of the phase transition from spin-stripe to spin-antiferro state 
as seen in Fig. [9] 



Let us try to understand the phase diagram of Fig. [9] 
in terms of these coupling constants. We begin with 
the small-77 regime where the groundstate is spin-stripe 
and orbital-ferro. Approaching from the high tempera- 
ture limit, we should look at the spin (orbital) couplings 
with disordered orbitals (spins). The strongest constants 
turn out to be Jd, which is slightly larger than J^ and 
Jy, all three being antiferromagnetic. This suggests that 
the system should undergo a transition to a spin-stripe 
state consistent with the phase diagram. The coupling 
constants of the orbital model are much weaker in the 
small-77 regime. The largest constant is K^^ suggesting 
an orbital-stripe order. However, since the spin- stripe 
state sets in at higher temperatures, in order to deter- 
mine the orbital order one should look at the coupling 
constants corresponding to the spin-stripe state. There 
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FIG. 11: (Color online) Effective exchange couplings Jx, Jy 
for nearest-neighbor and Jd for next-nearest-neighbor spins as 
a function of r? = Jh /U for 5 = 0.01, 7 = -0.05, and A = 0.7. 
The solid lines correspond to the couplings resulting for the 
corresponding orbital ground states whereas the dashed lines 
correspond to the orbitally disordered case. Note that for the 
orbitally disordered case Jx = Jy for all values of -q. The 
vertical dashed line indicates the location in 77 of the various 
phase transition as seen in Fig. |9l 



are three main effects (compare the solid and dashed lines 
in the low rj regime in Fig. llOp , (i) x- and y-directions be- 
come inequivalent in the sense that the couplings along 
X are suppressed while those along y are enhanced, (ii) 
the diagonal couplings are reduced strongly, and (iii) a 
single-site term is generated which acts as magnetic field 
for the orbital pseudospins. It is in fact this single site 
term that controls the ordering of the orbitals at low tem- 
peratures. This also explains the qualitatively different 
behavior of the orbital-ferro order parameter observed 
in Fig. [8^. Within the spin-stripe order, the single-site 
term becomes weaker with increasing 77 whereas the di- 
agonal term increases. This leads to a transition in the 
orbital sector from an orbital-ferro to an orbital-stripe 
phase near 77 = 0.11. The region between 0.14 and 0.2 
in eta is very interesting. Approaching from the high 
temperature the spins order into the 'stripe' state but 
as soon as the orbitals order into stripe state at lower 
temperature the diagonal couplings Jd are strongly re- 
duced and become smaller than Jy/2. This destabilizes 
the spin-stripe state and leads to an spin antiferro order- 
ing. For larger 77 the orbital ordering occurs at higher 
temperature. There is another transition slightly below 
rj = 0.3 where spins order into a ferro state. This is sim- 
ply understood as Jx = — Jy from the coupling constants 
of the Heisenberg model. 



IV. MAGNETIC EXCITATION SPECTRA 

We now set out to compute the magnetic excitation 
spectra, treating the orbital pseudospins as classical and 
static variables. Fixing the orbital degrees of freedom for 
a given set of parameters by the corresponding ground- 
state configuration we are left with an 5 = 1 Heisen- 
berg model written in (|18p . The exchange couplings are 
plotted in Fig. [TlJ Assuming the presence of local mo- 
ments, such J1-J2 models with a sufficiently large next- 
nearest neighbor exchange have been motivated and used 
to rationalize the (tt, 0) magnetism in the iron pnictidcs — 
and been used subsequently to calculate the magnetic 
excitation spectr a^^'^^ , where the incorporation of a rel- 
atively strong anisotropy between the nearest-neighbor 
couplings turned out to be necessary to understand the 
low energy spin- wave excitations^. 

In the presence of orbital ordering such an anisotropy 
of the effective magnetic exchange couplings appears nat- 
urally. Both, the orbital ferro and the orbital stripe order 
lead to a seizable anisotropy in the nearest-neighbor cou- 
plings Jx, Jy, where the anisotropy is much stronger for 
the orbital ferro order (see Fig. [TT|) . An even more drastic 
effect is the huge suppression of Jd in the orbital-stripe 
regime. 

On a classical level, the magnetic transitions are easily 
understood in the spin-only model as discussed be- 
fore. The transition from the stripe- AF to the (tt, tt) AF 
at ?7 w 0.14 occurs exactly at the point where Jy = 2Jd 
whereas the transition from (tt, tt) to ferromagnetic order 
at 77 « 0.3 corresponds to the point Jx = — Jy 

We proceed to calculate the magnetic excitation spec- 
tra in the Q = (tt, 0) and (tt, tt) phases within a linear 
spin-wave approximation. The classical ground states are 
given given by Sr = S{0,0,ar) with ar = exp(zQr) — 
±1. After performing a simple spin rotation, = 
— (JrS^, and = (7rS^, we express the ro- 
tated spin operators by Holstcin-Primakoff bosons, 5+ = 
V2S - hb, S~ = 5tV25- n, and 5^ = S-n with n = b^b 
to obtain the spin-wave Hamiltonian 

n = sl [Agib% + b^gbl^) + Bgiblbl^ + b^^bg)}, 



with 



Bq = J; 



T ^ 1 + COsQa; 
-JxCOsQx + Jx ^ COS Qx 

' Jd COS Qx COS Qy 

h-^(l + cos Qx COS Qy) COS Qx COS Qy 

1 — cos Qx , 1 — cos Qy 

COSQx + Jy C0S( 

J(i(l — COS Qx COS Qy) COS qx cos qy, 



X ^ y 



yielding the spin- wave dispersion Uq = S 
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FIG. 12: (Color online) Spin- wave excitation spectra for dif- 
ferent values of rj. Top: (tt, 0) magnet for disordered orbitals. 
The spectral weights are coded by line thickness and color, 
high intensity corresponds to red, low intensity to blue. Mid- 
dle: (7r,0) magnet for orbital-ferro {rj = 0.03,0.07,0.11) and 
orbital-stripe order (r; = 0.12,0.13,0.14). Bottom: (tt, tt) 
magnet with orbital-stripe order {ri = 0.15, 0.17, . . .). 

the inelastic structure factor at zero temperature^ 



Sinci{q,u;) = J—^S{u;-u;g) (19) 
V l+7q 

with 7q = Bq/Aq. Thc resulting excitation spectra 
are shown in Fig. [12] for different values of rj. In the 
case of disordered orbitals, the (tt, 0) antiferromagnet 
order is stable up to 77 « 0.25. Since Jx — Jy the 
spectrum is gapless not only at the ordering wave vec- 
tor (tt, 0) but also at the antiferromagnetic wave vector 
(tt, tt). However, the spectral weight is centered close 
to the ordering wave vector and goes strictly to zero 
at the antiferromagnetic wave vector. In the presence 
of orbital ordering the next-nearest neighbor couplings 
are anisotropic > Jy which in the case of the (tt, 0)- 
AF leads to a gap at the antiferromagnetic wave vector, 
tt) = 2-y(2Jc( — Jy){Jx — Jy)- Since the anisotropy 



and the diagonal exchange are large in the orbital ferro 
state we find a very big gap at (tt, tt). This gap reduces 
drastically for bigger 77 where the orbital stripe state be- 
comes favorable. Due to the large reduction of Jd and 
also of the anisotropy, the gap is considerably smaller 
and continuously goes to zero as we approach the tran- 
sition to the (tt, 7r)-AF at 77 « 0.14 where 2Jd — Jy = 0. 
This of course also leads to a strong anisotropy of the 
spin-wave velocities, Vy/vx = \/ {2Jd — Jy)/{2Jd + Jx)- 
On approaching the magnetic transition we find a signif- 
icant softening of modes along the (tt, 0) — (tt, tt) direction 
which leads to a considerable reduction of magnetic mo- 
ments close to the transition. 



V. DISCUSSION AND CONCLUSIONS 

In the preceding we have derived and studied a spin- 
orbital Kugel-Khomskii Hamiltonian relevant to the Fe- 
As planes of the parent compound of the iron supercon- 
ductors. A variety of interesting spin and orbital or- 
dered phases exist over a physical regime in parameter 
space. Due to the peculiarities of the pnictide lattice 
and this particular crystal field state we show that the 
relevant Kugel-Khomskii model is of a particularly inter- 
esting kind. 

The essence of the 'spin-charge-orbital' physics is dy- 
namical frustration. With so many 'wheels in the equa- 
tion' it tends to be difficult to find solutions that satisfy 
simultaneously the desires of the various types of degrees 
of freedom in the problem. This principle underlies the 
quite complex phase diagrams of for instance mangan- 
ites. But this dynamical frustration is also a generic 
property of the spin-orbital models describing the Jahn- 
Teller degenerate Mott-insulators. In the 'classic' Kugel- 
Khomskii model^S describing eg degenerate S" = 1/2 Sd^ 
systems of cubic 3d systems, Feiner et al^-^ discovered 
a point in parameter space where on the classical level 
this frustration becomes perfect. In the present context 
of pnictides this appears as particularly relevant since 
this opens up the possibility that quantum fluctuations 
can become quite important. 

We propose two specific orbital ordered phases that ex- 
plain the orthorhombic transition observed in the exper- 
iments. These are orbital-ferro and orbital-stripe states. 
The orbital-stripe order is particularly interesting since it 
leads to a spin model that provides possible explanation 
for the reduction of magnetic moment. It is our main 
finding that in the idealized pnictide spin-orbital model 
the conditions appear optimal for the frustration physics 
to take over. We find large areas in parameter space 
where frustration is near perfect. The cause turns out 
to be a mix of intrinsic frustration associated with hav- 
ing t2g type orbital degeneracy, and the frustration of a 
geometrical origin coming from the pnictide lattice with 
its competing " Jl — J2" superexchange pathways. The 
significance of this finding is that this generic frustration 
will render the spin-orbital degrees of freedom to be ex- 
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tremely soft, opening up the possibility that they turn 
into strongly fluctuating degrees of freedom ~ a desired 
property when one considers pnictide physics. 

We argued that the orthorhombic transition in half 
filled pnicitides and the associated anomalies in transport 
properties can be related to orbital order. When the pa- 
rameters are tuned away from the frustration regime the 
main tendency of the system is to anti-ferro orbital order- 
ing, which is the usual situation for antiferromagnets. An 
important result is that in the regime of relevance to the 
pnictides where the frustrations dominate we find phases 
that are at the same time (tt, 0) magnets and forms of or- 
bital order that are compatible with orthorhombic lattice 
distortions (Fig.'s 4,5). Besides the literal ferro-orbital 
ordered state (Fig. 5a), we find also a (tt, 0) or 'tweed' 
orbital order (Fig. 5b). This appears to be the more 
natural possibility in the insulating limit and if the weak 
superlattice reflections associated with this state would 
be observed this could be considered as a strong support 
for the literalness of the strong coupling limit. Surely, the 
effects of itinerancy are expected to modify the picture 
substantially. Propagating fermions are expected to sta- 
bilize ferro-orbital order o^^'"^^ , which enhances the spatial 
anisotropy of the spin-spin interactions further—. 

Among the observable consequences of this orbital 
physics is its impact of the spin fluctuations. We con- 
clude the paper with an analysis of the spin waves in 
the orbital ordered phases, coming to the conclusion that 
also the spin sector is quite frustrated, indicating that the 
quantum spin fluctuations should be quite strong offering 
a rational for a strong reduction of the order parameter. 

Thus we have forwarded the hypothesis that the un- 
doped iron pnictides are controlled by a very similar 
'spin-charge-orbital' physics as found in ruthenates and 
manganites. To develop a more quantitative theoretical 
expectation is less straightforward and as it is certainly 
beyond standard LDA and LDA-I-U approaches will re- 
quire investigations of correlated electron models such as 
we have derived here^i^, taking note of the fact that 
the pnictides most likely belong to the border line cases 
where the Hubbard U is neither small nor large compared 
to the bandwidth^. 



APPENDIX A: EFFECTIVE INTERACTION 
AMPLITUDES 



By acting with the hopping operator Tit (HI) on all 
states in the ground state sector and calculating the over- 
lap of the resulting states projected on the different in- 
termediate states wc find the effective interaction ampli- 
tudes. For the high-spin intermediate state (ri = 1) we 
find by projecting on the intermediate multiplet. 



f(i) 



J++ 
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J ZX 

f(i) 

J z 

f(l) 

J X 



(1) 



<b 
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^tab{tbb — taa) 



be 



e(4^2) + A 

^£(4^2) + A 

1 2tl, + tL + 1 



bb 



j-2 1 j-2 



be 



ei'A,) 



(Al) 
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where the hopping matrix elements have to be specified 
for a particular bond. Likewise, we find by projections 
on the intermediate low-spin states (n = 2) 
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The terms bilinear in the pseudospin operators resuh 
solely from hopping processes involving the \a) and \b) 
orbitals only. The hoppings between the |c) :— \xy) or- 
bitals enter only as a positive constant in Q'^2) leading 
to a conventional antiferromagnetic superexchange con- 
tribution. Interestingly, the coupling between the |c) and 
I a), 1 6) orbitals results in magnetic field terms for the or- 
bital pseudospins. 



x,yz 
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m2 + n2 = 1) with 



Using direction cosines (/, m, n) {P 
\l\ = I TO I resulting from the orthorhombic symmetry, we 
find that only the following hopping-matrix elements are 
non-zero, 
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These hopping matrix elements which are shown 
schematically in Fig.[2]can be parametrized by the lattice 
parameter A = \n/l\ and the ratio 7 = {pdT:)/(pda) as 



h/t 

t2/t 

hit 
hit 
hit 
hit 
hit 



-2(^2 - ^2 - C2) 

-{B^ + - C^) 
2AB - C2 
2^2 

2(S/A)2- A2 

2{AC + ABIX- B^IX), 



(B3) 



where we have introduced the overall energy scale t 
(pda)'^ I Apd and defined for abbreviation 



APPENDIX B: HOPPING MATRIX ELEMENTS 



For a given As-Fe bond with direction cosines l,Tn,n, 
the p to t2g hoppings are given by^ 



A 



B 



A(V3-27) 



C = 



V2 + y 
A(y3 + A27) 



V2TA^ 
^/3A2 + (2-A2)7 

V2TA^' 
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APPENDIX C: ORBITAL PART OF THE 
HAMILTONIAN 

For the nearest neighbor bonds along x and y the or- 
bital operators in the spin-orbital Hamiltonian are given 

by 



9\ 



92 



1 -h 2r? - (52 

1 + 7? + 2^ 
l + 2r}{l + 5) + 25' 
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-^a7(n + 7f2 + 3f3 + 35i)- ia^52, (C2) 

where we have defined fi = 'drj + d), f2 = fi\jj=oj 

rs = 1/(1 + 2?7 + (5), and rj = ri\s=o and introduced the 
functions 



Likewise, for the bonds along the x ±y diagonals we 
obtain 
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